Colossal magnetoresistance in manganites as a multicritical phenomenon 
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The colossal magnetoresistance in manganites j4Mn03 is studied from the viewpoint of mul- 
ticritical phenomena. To understand the complicated interplay of various phases, we study the 
Ginzburg-Landau theory in terms of both the mean-field approximation and the renormalization- 
group analysis to compare with the observed phase diagram. Several novel features, such as the 
first-order ferromagnetic transition, and the dip in the transition temperature near the multicritical 
point, can be understood as driven by enhanced fluctuations near the multicritical point. Further- 
more, we obtain a universal scaling relation for the H/M-M 2 plot (Arrott plot), which fits rather 
well with the experimental data, providing the further evidence for the enhanced fluctuation. 

PACS numbers: 75.30.Vn 75.40.Cx 64.60.Kw 75.60.Ej 
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Colossal magnetoresistance (CMR) in manganites is 
one of the most dramatic phenomena shown by strongly 
correlated electronic systems, and extensive experimental 
studies have revealed many aspects of this effect [[jj ||. 
However, its mechanism has been the subject of long- 
standing debates; many theories have been proposed such 
as double exchange j|, polaronic effect 0], and phase 
separation combined with percolation |5[ |(|and Griffiths 
singularity Shown in Figs. |l]are the phase diagrams 
of the CMR manganites. They clearly evidence that 
the CMR is related to the concomitant antiferromagnetic 
(AF) spin ordering, charge ordering (CO), and orbital or- 
dering (00) around the hole concentration x = 0.5. Near 
the phase boundary between the AF/CO/OO and the 
ferromagnetic metallic (FM) state, the transition tem- 
perature has a sharp dip and the critical magnetic field 
H c is reduced considerably. Hence, the CMR is collective 
in nature and differs from single particle properties such 
as the transition from small to large polarons. 

One of the subtle issues is the effect of the randomness 
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FIG. 1: Phase diagram for the mangan- 

ites: (i) Pro.55(Cai_j / Sr j ,)o.45Mn03 § and (ii) 
(Ndi_j / Sm !/ )o.55Sro.45Mn03 M. Solid squares represent 
a first-order transition. Red dots represent the data used 
in the scaling analysis in the text. The red region at 
OK < T < 50K in (ii) is a spin-glass phase. 



on the CMR. The ramdomness induced by the alloying 
and/or the nonuniform strain affects physical quantities 
near the boundary of the two phases, and alters the crit- 
ical properties easily. Figures |l| show the two systems 

(i) Pro. 55 (Cai_j / Sry)o.45Mn0 3 and (ii) (Ndi_j / Sm J/ )o.55 
Sro.45Mn03. It is not trivial which system is more dis- 
ordered, namely alloying (i) alkaline-earth atoms or (ii) 
rare-earth atoms. There are, however, three reasons to 
believe that the system (i) shows more intrinsic proper- 
ties than (ii) which is mostly dominated by the disorder 
effect. The first is that the strong suppression of the fer- 
romagnetic transition temperature and the appearance 
of the spin-glass type state in (ii) is well reproduced 
by recent model calculation in ||. The second is that 
the phase diagrams in ordered and disordered half-doped 
manganites Ln^Ba^MnOa [[l0| [IT] closely resemble 
Fig. [I] (i) and (ii), respectively. The last one is that the 
scaling fit works almost perfectly in (i), while it does only 
in the limited region in (ii), as we show in this Letter. In 

(ii) , the region near the phase boundary is dominated by 
the disorder effect. An appealing scenario for this behav- 
ior is the phase separation and/or percolation of the con- 
ducting paths H . This assumes a mixture of metallic and 
insulating domains; by applying magnetic field the metal- 
lic domains expand to result in the CMR. This is the 
static picture of the resistance network model controlled 
by the magnetic field. However, experiments on diffuse 
X-ray scattering and Raman scattering revealed that in- 
trinsic fluctuation is dynamic in nature Therefore, we 
need to take into account thermal and/or quantum fluc- 
tuations near the phase boundary. It is well-known that 
fluctuation is enhanced near a multicritical point, where 
more than two orders collide with each other. In this 
Letter we propose a new scenario based on this idea that 
the CMR originates from the enhanced fluctuation near 
the multicritical point, which is controlled by the small 
external magnetic field. We construct the phenomeno- 
logical Ginzburg-Landau (GL) model based on symmetry 
argument, and next give a renormalization-group (RG) 
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analysis for the multicritical phenomena to compare with 
experiments. This picture explains the scaling law and 
the first-order ferromagnetic phase transition as well as 
the enhanced sensitivity to the external magnetic field 
near the phase boundary. 

The ordering pattern of AF/CO/OO is complicated 
with an enlarged unit cell. Several microscopic mod- 
els for this state have been proposed Q. We employ 
here instead the phenomenological GL theory. We clas- 
sify possible terms in the free energy functional according 
to the symmetry of the order parameters. The relevant 
order parameters are those of ferromagnetism M, antifer- 
romagnetism S, charge ordering p, and orbital ordering 
T . Here, we discuss the dimensionality of each order pa- 
rameter. Both M and S are three-dimensional, while p 
is scalar. The orbital pseudo-vector T is originally three- 
dimensional but in the presence of the Jahn- Teller inter- 
action, which prefers real linear combinations of the two 
wavefunctions x 2 — y 2 and 3z 2 — r 2 , it should be regarded 
as two-dimensional: T = (T X ,T Z ). The free-energy func- 
tional should be rotationally invariant in the spin space, 
but not in the orbital pseudo-spin space. Hence, M and 
S should appear in the form of M 2 and S 2 . In contrast, 
third-order terms in T's are allowed, because —T is not 
iquivalent to — T. Linear terms in T's are generated by 
the coupling to the external uniaxial strain field. 

Next we consider the wavenumbers of the order pa- 
rameters. From the spatial pattern of the AF/CO/OO, 
we can easily see that wavenumber of each order 
is the following, p: (n, tt, 0), S: (tt,0, 7r), (0, 7r, 7r), 
(tt/2, -tt/2, tt), (-tt/2, tt/2, tt), f: (0,0,0), (tt,^0), 
(7r/2,7r/2,0), (-tt/2, -tt/2, 0), M: (0,0,0). Therefore, 
the allowed terms in the GL functional are 



F = 



(VM) 2 + (VS) 2 + (VT ) 2 + {Vpf 



+r M (M) 2 + r S (S) 2 + r T (f) 2 + r p (pf 

+ (r x pT T x + r z pT T z )p + g pS pS 2 + <? T Re(T x + iT z f 

+u M ((M) 2 ) 2 + u s ((S) 2 ) 2 + u T ({f) 2 ) 2 + u pP 4 



(1) 



Here some remarks are in order. First, the bilinear cou- 
pling between p and T enforces that these two orders 
accompany each other. This agrees with experiments. 
Second, the third-order term in T make the transition 
first-order. However the magnitude of the jump at the 
transition depends on the relative values of r p and ry; 
when r p <C rx, it becomes nearly second-order, while it 
is strongly first-order in the other limit. Experimentally, 
the CO/OO transition in the narrow band- width side is 
nearly second-order, which means that the transition is 
mainly driven by the CO. Hence, we neglect below the 
orbital ordering T. Third, the term pS 2 prohibits the AF 
without the CO, namely the CO occurs at a temperature 
higher than or as high as the AF. This is a generic con- 
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FIG. 2: Mean- field phase diagrams for the GL functional (|f|) 
for in the W-T plane. Here W is the bandwidth and T is the 
temperature. The broken (solid) lines represent the second 
(first)-order phase transitions. 



elusion from the symmetry independent of microscopic 
mechanisms, and is consistent with the experiments. 

Figure || shows three possible mean-field phase dia- 
grams for the GL functional ([!]), in the plane of the 
band-width W and the temperature T. Figure ||(a) is 
the most relevant to the experimental results. Nonethe- 
less, this mean-field analysis cannot capture several ex- 
perimental features. One is the first-order FM transi- 
tion in the wide-bandwidth side. Another is the dip of 
the transition temperatures of FM and CO/OO near the 
critical bandwidth W c . These two features are due to 
the fluctuations enhanced near W c . This is the so-called 
multicritical phenomenon between the competing orders. 

We now turn to the RG analysis of this fluctuation. 
Let us consider a system with competing two orders "A" 
and "B". We express their order parameters as an N^- 
dimensional vector (\>a and an TVe-dimensional vector 
Assuming a rotational invariance in §a and in we can 
write its classical free energy up to quartic order as 

F cl = ^{r A <t> 2 A + r B cf> 2 B ) + ^{udp A + 2w<\> 2 A <\> 2 B + v<f>%). 

(2) 

In the present case, the phases A and B correspond to the 
FM phase and the CO/OO phase. Thus, <j) A is a mag- 
netization M. The average value Na ^ a a ^ b J b is roughly 
proportional to the temperature measured from the mul- 
ticritical point, while r A — r B is proportional to the de- 
viation of bandwidth from the multicritical point. The 
model (||) was studied in Jll| by the RG. The fluctuation- 
induced first-order transition fl4|| is known to occur in 
(Q) when w > y/uv fll5|| . For an understanding of the 
fluctuation-induced first-order transition, an effective po- 
tential [|l6) is a useful tool. The effective potential is 
given by adding fluctuation contribution to the classical 
Lagrangian. Bare coefficients in the effective potential 
should be renormalized to remove ultraviolet divergence 
at d = 4. By adding counterterms as in jl6|], we get 

F = \{r A <t>\ + r B 4> 2 B ) + ^(wfa + 2w^ 2 B + v<j>%) 
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where 

7+ + 7- - + 7+7- = ftt/fe - 4to$^6) 

and k is a parameter setting a scale of momentum. We 
neglected the terms higher than 0(4 — d). The quantities 
7"i, «, i>, it; in (^) are renormalized ones and are finite. 

Let us consider the fluctuation-induced first-order 
transition [[[i] between the disordered phase and the "A" 
phase. We thus assume that the "B" field 4>b is not or- 
dered, i.e. 4>b — 0- In the RG language, the fluctuation- 
induced first-order transition occurs when the RG flow 
runs into an unstable region of the model (0). ft means 
that sixth-order terms, though omitted in (0), are nec- 
essary for stability, leading to a first-order transition. 
Hence, if the RG flow crosses the boundary of the stabil- 
ity region, u — 0, the system undergoes a first-order tran- 
sition to the "A" phase. In other words, we follow the RG 
flow for u, v, w |13[ ] until we reach the line u = 0. Let K\ 
denote the value of k when u(k) = 0. Other quantities 
as r a,tb,^>a are also renormalized. It is worth noting 
that these quantities are renormalized multiplicatively, 
i.e. n(Ki) = nf(Ki), (/)a{ki) = cj>Ag{n\), where n and 
4>a are the initial values. Thus, renormalization is merely 
a change of scale for them. The condition m(«i) = 
greatly simplifies the above free energy as 

F = \ta<\>\ + N A f(r A ) + N B f(r B + w<j? A ), (7) 

where f(x) = ^- (In a; + |) and Ki is set as unity since it 
can be absorbed by the change of scale of other variables. 
The equation of state is 

dF 

H A = tt- = r A (pA + 24> A N B wf{r B +W<f> 2 A ). (8) 

0(f>A 

where Ha is a field conjugate to the 4>A field. Since we 
identify this 4>A as a magnetization M, Ha is a magnetic 
field H. It is convenient to rewrite (ft) as 

^ = t a + 2wN B f(r B + wM 2 ). (9) 
M 

Therefore if we plot H/M versus M 2 at various temper- 
ature and bandwidth, all the plots will be degenerate. In 
other words, the H/M-M 2 curve, called the Arrott plot, 
will undergo parallel transport when we change temper- 
ature or bandwidth. One would notice that in a cer- 
tain range of temperature and bandwidth this universal 
curve crosses the horizontal axis once or twice. In that 
case some part of the universal curve becomes unphysical, 
and the system undergoes a first-order phase transition 
by increasing the magnetic field. The system is either 
ferromagnetic or metamagnetic. 



To verify this scenario, we used two series of data in 
Figs. (i) for y = 0.25, 0.3, 0.4 and (ii) for 0.1 < y < 0.4, 
with varying y and T. We can regard y as a parameter 
controlling the bandwidth. We expect that ta and tb are 
functions of y and T while w is a constant. We expand 
ta and in the vicinity of the multicritical point as 
TA = cat AT + CAy&y, r B = c BT AT + c B yAy, where 
cat, cbt, CAy, CB y are constants. In view of (||), cat, cbt 
(cAy, CBy) represents an amount of parallel displacement 
of the H/M-M 2 plot when the temperature T (doping y) 
is changed. We fitted the data in the following way. First 
we varied cat, cbt, CAy, cb v so that the plots for various y 
and T overlap most after the parallel displacement. Then 
the scales of the abscissa and the ordinate are varied to fit 
to (H). We used the data shown as red dots in Figs. [|. We 
discarded the data with M > 1.6/ib, because when the 
magnetization approaches the saturation value, the GL 
functional up to quadratic order is no longer appropriate. 
We also discarded the data for 0.5 < y < 0.8 in (ii). They 
do not fit well with the scaling curve. This is reasonable 
because in (ii) critical fluctuations are washed out near 
the multicritical point, by which the scaling plot deviates 
from (||). 

The resulting plot is shown in Fig. S. The plots fall 
into one curve with good accuracy. This shows that the 
fluctuation is enhanced near the multicritical point. Be- 
cause the critical scaling @ holds for a wide range of 
data shown in Fig. ||, the critical region is rather large 
(~ 80K). It manifests strong correlation of electrons in 
the manganites. Whether the system is critical or not is 
determined by £, measured by a length scale £q ~ a , 
where A is a gap, a is a lattice constant, and E-p is 
the Fermi energy. In the manganites A ~ E-p yields 
£o ~ o,. Temperature dependence of £ in a related com- 
pound Pr .5Cao.5Mn03 can be obtained from Fig. 2(c) 
of |0. It is i ~ 100A, near the transition temperature 
T c = 235K, and is £ ~ 20A even at 300K, which is 65K 
higher than T c . Therefore, £ is much longer that £o in 
the wide temperature range, implying that the critical 
region should be large. 

The enhanced fluctuation makes the system sensitive 
to the external magnetic field H. The CO/OO state 
easily becomes the FM state by a weak magnetic field. 
Let H c denote this critical field. An exponent x defined 
by H c (x (W c — W) x represents the sensitivity to the 
external field. Without the fluctuation, x should be equal 
to unity. On the other hand, in the multicritical region 
but not in the fluctuation-induced first-order-transition 
region, x is controlled by the bicritical fixed point, and 
is larger than unity. Calculation based on the (4 — e)- 
expansion up to 0(e 2 ) results in x — ~ 1.37, where <fi 
is the crossover exponent. It can be compared with the 
experimental value x ~ 1.6, extracted from the data for 
Pro.65(Cai_j / Sr J/ )o.35Mn03 Q. This increase of x near 
the multicritical point, i.e. sensitivity to the external 
field, emerges as the CMR phenomena. 
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temperature limit, dynamics and statics are coupled and 
the Griffiths singularities influence the whole quantum 
critical phenomena. In the manganites, this possibility 
seems to be prevented by the glassy state appearing at 
low temperatures. 
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FIG. 3: Scaling plot for (i) Pro.55(Cai_j / Sr B )o.45Mn03 and 
(ii) (Ndi_ !/ Smj / )o.55Sro.45Mn03. The broken line represent 
the f'(x) — j (lnx + 1) plot. Each dot show the data for fixed 
T and y. 

Here we discuss the effect of disorder, and the rele- 
vance of the Griffiths singularity j?]] to the present study. 
It has been pointed out that in the temperature region 
T c < T < Tc°^ where T^ ' is the transition temperature 
for the pure system while T c is the suppressed one by 
the random dilution, there appears some singularity in 
the time-dependence of correlation functions of the or- 
der parameter, such as the extended power-law decay. 
However, this singularity does not occur in the thermo- 
dynamic quantities which we discussed in this Letter. It 
is possible that the fixed point for the pure system is un- 
stable against disorder, but according to the Harris crite- 
rion the bicritical fixed point described above is at least 
locally stable . Hence, the thermodynamic properties 
are described by the usual bicritical fixed point, although 
the location of the critical point is shifted. However, the 
disorder effect on the first-order phase transition are be- 
yond the scope of the present work, and ref. Q is rele- 
vant to this issue. The percolation scenario could also 
lead to the CMR, where the resistivity will depend sen- 
sitively on random realization of the metallic paths etc., 
and consequently on samples. This behavior has been 
actually observed in dilutely Cr-doped manganites [fl9|| . 
Cr ions destroy the CO/OO locally, and introduce the 
FM region. In these dilutely doped samples, the resistiv- 
ity at low temperature depends on the Cr-concentration 
dramatically, and the hysteresis appears in the temper- 
ature cycle. In contrast, when the Cr-concentration in- 
creases, the resisivity no longer depends on samples or 
heat cycle. In the latter case, the thermodynamic phases 
are well-defined and the CMR is triggered by the phase 
change between them, which is the subject of our study 
here. Hence, there are two types of the CMR; one is due 
to the percolating path and the other is due to the mul- 
ticritical fluctuation near the phase change. In the zero 
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